Utility of polygenic risk scores in UK cancer screening: a modelling analysis.

BACKGROUND: It is proposed that, through restriction to individuals delineated as high risk, polygenic risk scores (PRSs) might enable more efficient targeting of existing cancer screening programmes and enable extension into new age ranges and disease types. To address this proposition, we present an overview of the performance of PRS tools (ie, models and sets of single nucleotide polymorphisms) alongside harms and benefits of PRS-stratified cancer screening for eight example cancers (breast, prostate, colorectal, pancreas, ovary, kidney, lung, and testicular cancer). METHODS: For this modelling analysis, we used age-stratified cancer incidences for the UK population from the National Cancer Registration Dataset (2016-18) and published estimates of the area under the receiver operating characteristic curve for current, future, and optimised PRS for each of the eight cancer types. For each of five PRS-defined high-risk quantiles (ie, the top 50%, 20%, 10%, 5%, and 1%) and according to each of the three PRS tools (ie, current, future, and optimised) for the eight cancers, we calculated the relative proportion of cancers arising, the odds ratios of a cancer arising compared with the UK population average, and the lifetime cancer risk. We examined maximal attainable rates of cancer detection by age stratum from combining PRS-based stratification with cancer screening tools and modelled the maximal impact on cancer-specific survival of hypothetical new UK programmes of PRS-stratified screening. FINDINGS: The PRS-defined high-risk quintile (20%) of the population was estimated to capture 37% of breast cancer cases, 46% of prostate cancer cases, 34% of colorectal cancer cases, 29% of pancreatic cancer cases, 26% of ovarian cancer cases, 22% of renal cancer cases, 26% of lung cancer cases, and 47% of testicular cancer cases. Extending UK screening programmes to a PRS-defined high-risk quintile including people aged 40-49 years for breast cancer, 50-59 years for colorectal cancer, and 60-69 years for prostate cancer has the potential to avert, respectively, a maximum of 102, 188, and 158 deaths annually. Unstratified screening of the full population aged 48-49 years for breast cancer, 58-59 years for colorectal cancer, and 68-69 years for prostate cancer would use equivalent resources and avert, respectively, an estimated maximum of 80, 155, and 95 deaths annually. These maximal modelled numbers will be substantially attenuated by incomplete population uptake of PRS profiling and cancer screening, interval cancers, non-European ancestry, and other factors. INTERPRETATION: Under favourable assumptions, our modelling suggests modest potential efficiency gain in cancer case detection and deaths averted for hypothetical new PRS-stratified screening programmes for breast, prostate, and colorectal cancer. Restriction of screening to high-risk quantiles means many or most incident cancers will arise in those assigned as being low-risk. To quantify real-world clinical impact, costs, and harms, UK-specific cluster-randomised trials are required. FUNDING: The Wellcome Trust.

Characterization of a mutant samhd1 zebrafish model implicates dysregulation of cholesterol biosynthesis in Aicardi-Goutières syndrome.

Aicardi-Goutières syndrome (AGS1-9) is a genetically determined encephalopathy that falls under the type I interferonopathy disease class, characterized by excessive type I interferon (IFN-I) activity, coupled with upregulation of IFN-stimulated genes (ISGs), which can be explained by the vital role these proteins play in self-non-self-discrimination. To date, few mouse models fully replicate the vast clinical phenotypes observed in AGS patients. Therefore, we investigated the use of zebrafish as an alternative species for generating a clinically relevant model of AGS. Using CRISPR-cas9 technology, we generated a stable mutant zebrafish line recapitulating AGS5, which arises from recessive mutations in SAMHD1. The resulting homozygous mutant zebrafish larvae possess a number of neurological phenotypes, exemplified by variable, but increased expression of several ISGs in the head region, a significant increase in brain cell death, microcephaly and locomotion deficits. A link between IFN-I signaling and cholesterol biosynthesis has been highlighted by others, but not previously implicated in the type I interferonopathies. Through assessment of neurovascular integrity and qPCR analysis we identified a significant dysregulation of cholesterol biosynthesis in the zebrafish model. Furthermore, dysregulation of cholesterol biosynthesis gene expression was also observed through RNA sequencing analysis of AGS patient whole blood. From this novel finding, we hypothesize that cholesterol dysregulation may play a role in AGS disease pathophysiology. Further experimentation will lend critical insight into the molecular pathophysiology of AGS and the potential links involving aberrant type I IFN signaling and cholesterol dysregulation.

Genome-wide association analysis identifies a susceptibility locus for sporadic vestibular schwannoma at 9p21.

Vestibular schwannomas are benign nerve sheath tumours that arise on the vestibulocochlear nerves. Vestibular schwannomas are known to occur in the context of tumour predisposition syndromes NF2-related and LZTR1-related schwannomatosis. However, the majority of vestibular schwannomas present sporadically without identification of germline pathogenic variants. To identify novel genetic associations with risk of vestibular schwannoma development, we conducted a genome-wide association study in a cohort of 911 sporadic vestibular schwannoma cases collated from the neurofibromatosis type 2 genetic testing service in the north-west of England, UK and 5500 control samples from the UK Biobank resource. One risk locus reached genome-wide significance in our association analysis (9p21.3, rs1556516, P = 1.47 × 10-13, odds ratio = 0.67, allele frequency = 0.52). 9p21.3 is a genome-wide association study association hotspot, and a number of genes are localized to this region, notably CDKN2B-AS1 and CDKN2A/B, also referred to as the INK4 locus. Dysregulation of gene products within the INK4 locus have been associated with multiple pathologies and the genes in this region have been observed to directly impact the expression of one another. Recurrent associations of the INK4 locus with components of well-described oncogenic pathways provides compelling evidence that the 9p21.3 region is truly associated with risk of vestibular schwannoma tumorigenesis.

Re-evaluation of missense variant classifications in NF2.

Missense variants in the NF2 gene result in variable NF2 disease presentation. Clinical classification of missense variants often represents a challenge, due to lack of evidence for pathogenicity and function. This study provides a summary of NF2 missense variants, with variant classifications based on currently available evidence. NF2 missense variants were collated from pathology-associated databases and existing literature. Association for Clinical Genomic Sciences Best Practice Guidelines (2020) were followed in the application of evidence for variant interpretation and classification. The majority of NF2 missense variants remain classified as variants of uncertain significance. However, NF2 missense variants identified in gnomAD occurred at a consistent rate across the gene, while variants compiled from pathology-associated databases displayed differing rates of variation by exon of NF2. The highest rate of NF2 disease-associated variants was observed in exon 7, while lower rates were observed toward the C-terminus of the NF2 protein, merlin. Further phenotypic information associated with variants, alongside variant-specific functional analysis, is necessary for more definitive variant interpretation. Our data identified differences in frequency of NF2 missense variants by exon between gnomAD population data and NF2 disease-associated variants, suggesting a potential genotype-phenotype correlation; further work is necessary to substantiate this.

MRSD: A quantitative approach for assessing suitability of RNA-seq in the investigation of mis-splicing in Mendelian disease.

Variable levels of gene expression between tissues complicates the use of RNA sequencing of patient biosamples to delineate the impact of genomic variants. Here, we describe a gene- and tissue-specific metric to inform the feasibility of RNA sequencing. This overcomes limitations of using expression values alone as a metric to predict RNA-sequencing utility. We have derived a metric, minimum required sequencing depth (MRSD), that estimates the depth of sequencing required from RNA sequencing to achieve user-specified sequencing coverage of a gene, transcript, or group of genes. We applied MRSD across four human biosamples: whole blood, lymphoblastoid cell lines (LCLs), skeletal muscle, and cultured fibroblasts. MRSD has high precision (90.1%-98.2%) and overcomes transcript region-specific sequencing biases. Applying MRSD scoring to established disease gene panels shows that fibroblasts, of these four biosamples, are the optimum source of RNA for 63.1% of gene panels. Using this approach, up to 67.8% of the variants of uncertain significance in ClinVar that are predicted to impact splicing could be assayed by RNA sequencing in at least one of the biosamples. We demonstrate the utility and benefits of MRSD as a metric to inform functional assessment of splicing aberrations, in particular in the context of Mendelian genetic disorders to improve diagnostic yield.

BBS1 branchpoint variant is associated with non-syndromic retinitis pigmentosa.

BACKGROUND: Inherited retinal diseases (IRDs) can be caused by variants in >270 genes. The Bardet-Biedl syndrome 1 (BBS1) gene is one of these genes and may be associated with syndromic and non-syndromic autosomal recessive retinitis pigmentosa (RP). Here, we identified a branchpoint variant in BBS1 and assessed its pathogenicity by in vitro functional analysis. METHODS: Whole genome sequencing was performed for three unrelated monoallelic BBS1 cases with non-syndromic RP. A fourth case received MGCM 105 gene panel analysis. Functional analysis using a midigene splice assay was performed for the putative pathogenic branchpoint variant in BBS1. After confirmation of its pathogenicity, patients were clinically re-evaluated, including assessment of non-ocular features of Bardet-Biedl syndrome. RESULTS: Clinical assessments of probands showed that all individuals displayed non-syndromic RP with macular involvement. Through detailed variant analysis and prioritisation, two pathogenic variants in BBS1, the most common missense variant, c.1169T>G (p.(Met390Arg)), and a branchpoint variant, c.592-21A>T, were identified. Segregation analysis confirmed that in all families, probands were compound heterozygous for c.1169T>G and c.592-21A>T. Functional analysis of the branchpoint variant revealed a complex splicing defect including exon 8 and exon 7/8 skipping, and partial in-frame deletion of exon 8. CONCLUSION: A putative severe branchpoint variant in BBS1, together with a mild missense variant, underlies non-syndromic RP in four unrelated individuals. To our knowledge, this is the first report of a pathogenic branchpoint variant in IRDs that results in a complex splice defect. In addition, this research highlights the importance of the analysis of non-coding regions in order to provide a conclusive molecular diagnosis.

Comparison of in silico strategies to prioritize rare genomic variants impacting RNA splicing for the diagnosis of genomic disorders.

The development of computational methods to assess pathogenicity of pre-messenger RNA splicing variants is critical for diagnosis of human disease. We assessed the capability of eight algorithms, and a consensus approach, to prioritize 249 variants of uncertain significance (VUSs) that underwent splicing functional analyses. The capability of algorithms to differentiate VUSs away from the immediate splice site as being 'pathogenic' or 'benign' is likely to have substantial impact on diagnostic testing. We show that SpliceAI is the best single strategy in this regard, but that combined usage of tools using a weighted approach can increase accuracy further. We incorporated prioritization strategies alongside diagnostic testing for rare disorders. We show that 15% of 2783 referred individuals carry rare variants expected to impact splicing that were not initially identified as 'pathogenic' or 'likely pathogenic'; one in five of these cases could lead to new or refined diagnoses.

Modelling the developmental spliceosomal craniofacial disorder Burn-McKeown syndrome using induced pluripotent stem cells.

The craniofacial developmental disorder Burn-McKeown Syndrome (BMKS) is caused by biallelic variants in the pre-messenger RNA splicing factor gene TXNL4A/DIB1. The majority of affected individuals with BMKS have a 34 base pair deletion in the promoter region of one allele of TXNL4A combined with a loss-of-function variant on the other allele, resulting in reduced TXNL4A expression. However, it is unclear how reduced expression of this ubiquitously expressed spliceosome protein results in craniofacial defects during development. Here we reprogrammed peripheral mononuclear blood cells from a BMKS patient and her unaffected mother into induced pluripotent stem cells (iPSCs) and differentiated the iPSCs into induced neural crest cells (iNCCs), the key cell type required for correct craniofacial development. BMKS patient-derived iPSCs proliferated more slowly than both mother- and unrelated control-derived iPSCs, and RNA-Seq analysis revealed significant differences in gene expression and alternative splicing. Patient iPSCs displayed defective differentiation into iNCCs compared to maternal and unrelated control iPSCs, in particular a delay in undergoing an epithelial-to-mesenchymal transition (EMT). RNA-Seq analysis of differentiated iNCCs revealed widespread gene expression changes and mis-splicing in genes relevant to craniofacial and embryonic development that highlight a dampened response to WNT signalling, the key pathway activated during iNCC differentiation. Furthermore, we identified the mis-splicing of TCF7L2 exon 4, a key gene in the WNT pathway, as a potential cause of the downregulated WNT response in patient cells. Additionally, mis-spliced genes shared common sequence properties such as length, branch point to 3' splice site (BPS-3'SS) distance and splice site strengths, suggesting that splicing of particular subsets of genes is particularly sensitive to changes in TXNL4A expression. Together, these data provide the first insight into how reduced TXNL4A expression in BMKS patients might compromise splicing and NCC function, resulting in defective craniofacial development in the embryo.

Machine Learning Approaches for the Prioritization of Genomic Variants Impacting Pre-mRNA Splicing.

Defects in pre-mRNA splicing are frequently a cause of Mendelian disease. Despite the advent of next-generation sequencing, allowing a deeper insight into a patient's variant landscape, the ability to characterize variants causing splicing defects has not progressed with the same speed. To address this, recent years have seen a sharp spike in the number of splice prediction tools leveraging machine learning approaches, leaving clinical geneticists with a plethora of choices for in silico analysis. In this review, some basic principles of machine learning are introduced in the context of genomics and splicing analysis. A critical comparative approach is then used to describe seven recent machine learning-based splice prediction tools, revealing highly diverse approaches and common caveats. We find that, although great progress has been made in producing specific and sensitive tools, there is still much scope for personalized approaches to prediction of variant impact on splicing. Such approaches may increase diagnostic yields and underpin improvements to patient care.

Disease modeling of core pre-mRNA splicing factor haploinsufficiency.

The craniofacial disorder mandibulofacial dysostosis Guion-Almeida type is caused by haploinsufficiency of the U5 snRNP gene EFTUD2/SNU114. However, it is unclear how reduced expression of this core pre-mRNA splicing factor leads to craniofacial defects. Here we use a CRISPR-Cas9 nickase strategy to generate a human EFTUD2-knockdown cell line and show that reduced expression of EFTUD2 leads to diminished proliferative ability of these cells, increased sensitivity to endoplasmic reticulum (ER) stress and the mis-expression of several genes involved in the ER stress response. RNA-Seq analysis of the EFTUD2-knockdown cell line revealed transcriptome-wide changes in gene expression, with an enrichment for genes associated with processes involved in craniofacial development. Additionally, our RNA-Seq data identified widespread mis-splicing in EFTUD2-knockdown cells. Analysis of the functional and physical characteristics of mis-spliced pre-mRNAs highlighted conserved properties, including length and splice site strengths, of retained introns and skipped exons in our disease model. We also identified enriched processes associated with the affected genes, including cell death, cell and organ morphology and embryonic development. Together, these data support a model in which EFTUD2 haploinsufficiency leads to the mis-splicing of a distinct subset of pre-mRNAs with a widespread effect on gene expression, including altering the expression of ER stress response genes and genes involved in the development of the craniofacial region. The increased burden of unfolded proteins in the ER resulting from mis-splicing would exceed the capacity of the defective ER stress response, inducing apoptosis in cranial neural crest cells that would result in craniofacial abnormalities during development.